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Abstract 

We study the 7-Z interference in the process e + e~ — > T — ► t + t~ 
if"} , as a means to measure the neutral current coupling of the b-quark. 

The helicity amplitudes are calculated from resonant and background 
diagrams and the spin density matrix of the final state is discussed. 
The spin analyzer of the r's is illustrated with the decays irv and 
^ \ pv — > (7T7t)za With 10 8 T a sensitivity to g v of a few per cent could be 

■ reachable. 
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1 Introduction 



A precise determination of the fermionic electroweak couplings can provide 
stronger hints on the nature of new physics at higher scales through the quan- 
tum corrections to the effective theory at lower energies. At present there 
is no solid experimental discrepancy with the predictions of the Standard 
Model 0, and this fact is used to set bounds on these new dynamics beyond 
the Standard Model. 

In this context, the agreement between the measured Z—ff couplings and 
the Standard Model is very good. However, while the lepton couplings have 
been independently measured for the three families at LEP, for the quarks 
only the b and c events can be separated in the hadronic event sample and 
consequently their couplings measured exclusively. From the measurements 
of Rb and Af, [0] we can get the values of the vector and axial couplings of 
Z to 66 |§] . To find the accuracy on the bottom-Z couplings obtainable from 
these measurements it is enough to use tree level expressions for Rb and Ab 
because higher order corrections will slightly shift the central values but will 
not modify the errors. Using the measurements quoted in 0, and defining the 
vector and axial couplings in the Standard Model as g\r = — 1/4 + sin 2 9 w /3 
and g\ = 1/4, the accuracy on the bottom-Z couplings is: Sgy = ±0.013 and 
5g A = ±0.007. It is important to notice that while the uncertainty in the g\ 
coupling is of 2.8%, the g\ measurent is worse, with an uncertainty of 7.5%. 

For the light quarks this separation has not been achieved at LEP and so 
an exclusive determination of their couplings is not available. A study of the 
final state distributions at different meson facilities can provide independent 
measurements of these couplings. In a previous work |4| , we showed that a $ 
factory with polarized e~ beams can supply the information on the s-quark 
couplings. 

We will show in this work, that a detailed study of the final state dis- 
tributions of the decay products of the r's for e + e~ — > T — > t + t~ would 
provide valuable information on the vector Z — bb coupling, g\. 

To determine this coupling at energies well below the Z pole, we will 
have to measure the interference between 7 and Z, where Z is coupled to a 
bb resonance. The bb mesons which can couple both to 7 and Z are the T 
mesons. We are interested in a process in which the T is coupled to a Z 
either in its production or in its decay. This means we will study the decay 
of a polarized T or its weak decay. 
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For these reasons, we will analyze the leptonic decays of the T resonances, 
therefore we could use T(1S), T(2S) and T(3S), but not T(45'), that decays 
dominantly to BB where the information on the T polarization is lost. In 
this work we will concentrate on the T(IS') resonance, but everything would 
be similar for T(2S) and T(3S). We can see in Table 1 that the branching 
ratios of T(IS') to the three charged leptons are approximately 3%, but e + e~ 
and yU + /i~ can not be used because their polarizations are not measurable in 
the detector through decay distributions. Obviously e + e~ are stable particles 
and [i + at this energy do not decay inside the detector. As we will see 
explicitely later, all the relevant information on g v , that comes from a P- 
odd 7 — Z interference, only appears at leading order in the polarization of 
the final leptons. This means that we are constrained to consider the decay 
T — > t + t~ and to measure the r polarizations. 

Before entering a complete analysis of the r observables, to estimate the 
sensitivity of this process to the vector coupling, let us first consider the r~ 
longitudinal polarization, suggested in Ref. || in another context. We will 
make a reasonable approximation in order to get a simple and clean result: 
the resonant diagrams will dominate the process on the T(1S) peak. So, we 
consider only diagrams (2), (3) and (4) in Fig. 1. Under these conditions we 
get a longitudinal polarization from the parity violating interferences of the 
dominant amplitude (2) with the neutral current amplitudes (3) + (4). 

P — ^ ( 1 + cos2e )\p\P° + 2cose (P°) 2 m 

z ' V2^a 9A Q b (l + cos 2 #)(p°) 2 + sin 2 #M2 U 

where g& = 1/4 is the axial coupling to the leptons, Qb the charge of the b 
quark and p^ = {p°,p{@)) the four-momentum of the r~. 

In this expression we can discover some interesting features of this ob- 
servable, 

• It is linear on g v , the vector coupling which we want to determine, 
showing up together with the axial coupling of the Z to leptons. 

• The magnitude of the r-polarization on the T peak is set by the factor 
8G F M|/(47ra v / 2) ^ 0.064, that translates into P z ,{9 = 0) = 0.032. 

• It is independent of the hadronic structure of the resonance which can- 
cels completely in the ratio. 
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All these properties will be modified when we include the non-resonant 
diagrams, but these new contributions will correct this result at the level of a 
few per cent, so this will be the dominant behaviour of our observables. These 
properties were already pointed out in Ref. ||, where it was shown that this 
polarization is enhanced in the vicinity of the resonance. Approximately at 
four amplitudes below the resonance one gets a polarization five times bigger, 
but the number of events decreases three orders of magnitude and so it is 
more efficient from the point of view of statistics to stay on the peak of the 
resonance. 

The above example indicates that a more detailed analysis of the problem 
is worth. In the next section we calculate the complete t~t + density matrix 
in this process which contains all the relevant information on g\r. In section 3 
we study the hadronic decays of the r to meassure the r density matrix, and 
from here we analyse, in section 4, the statistical accuracy that is possible to 
get in the meassure of the vector coupling to the b quark in each channel. 

2 t + t~ density matrix 

The density matrix from e~(Z_, £_)e + (/ + , £ + ) — > r~(jo_, A_)r + (jo + , A+), in 
terms of helicity amplitudes 0, and when our initial beams are unpolarized, 
Eq. ( |A.1| ), is given by, 

P(A_,A+),(A'_,A' + ) = Yl /(A_,A + ),r ( e )f(X'_,X' + ),C W ^ 
f 

where the angle (8) is given by the direction of the r~ relative to the initial 

e~ beam, with the x-z plane defined as the scattering plane. 

Using reduced helicity amplitudes, Ti , and rotation matrices these he- 
Ac; 

licity amplitudes are |7|], 

fuW = d U^ T k ( 3 ) 

where A = A_ — A + , £ = — £ + and dj x (9) are the reduced rotation matrices 
about the y-axis. If we neglect the electron mass, the total angular momen- 
tum of the process is always J = 1. Therefore we get Eq. (|3]) where we have 
re-defined our reduced helicity amplitudes with respect to Ref. 0, including 
in our definition several normalization factors irrelevant in our analysis ||]. 
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Furthermore, helicity conservation in the electron vertex implies that £ fixes 
completely £ = (£+,£_). These reduced helicity amplitudes get contribu- 
tions from diagrams (1) to (5) in Figure 1. For the dominant amplitudes 
contributing through interferences to P-odd and C-odd observables we can 
write 



T U = T a, 5 (7) + T u (<yZ A ) + T U {Z A1 ) + T U (Z AA ) (4) 

where T x ^(7) accounts for the contribution of diagrams (1) and (2). T x ^fZ A ) 
is the P-violating piece of diagram (3) plus the V A piece of diagram (5) and 
T^^(Z'-fA) the corresponding P-violating piece of diagram (4) plus the AV 
piece of diagram (5). Finally T^^(Za,a) is the contribution of diagram (5) 
with axial couplings in the initial and final vertices. Notice that here we have 
not included the T^^Zy), T^^(Zjv) and T^^Zyy) pieces, because they are 
sub-dominant with respect to T% ^(7), both in the resonant and non-resonant 
components. 

Taking into account the transformation properties under P of these am- 
plitudes we get 



T U W= T_ x> _ € ( 7 ) = 7^(7) =T_ X>€ (7) (5) 

T x,^ Z a) = ~T_i_^Z A ) = T x ^(jZ A ) = -T_ X ^Z A ) 

T^(Zai) = -T_ x ^(Zai) = -T^Zai) = T_^{Z A l) 

T^(Z A a) = T^^Zaa) = -T x _^Zaa) = -T_^(Z AA ) 

where the normalization of these reduced helicity amplitudes is such that 

^ = Tr(p°«') = 2sin 2 £|T ( ^ 

+ 4cos^e{r (+i _ )il ( 7 )r ( * +! _ )!l (Z^)} (6) 



(m + , +) ,i(7)i 2 +iT (+ ,-),i(7)r) (7) 

Notice that a C-odd forward-backward interference is generated by the axial 
couplings of the Z in both the electron and r vertices. 

We can calculate these reduced helicity amplitudes from the Feynman 
diagrams 1-5 of Fig. 1 following the method explained in App. B. 
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(8) 



where K is only a constant which takes care of the different normalization of 
the helicity amplitudes and the Feynman amplitudes. = — § is the charge 
of the b quark, gv(A) is the vector (axial) coupling to the leptons, P-f stands 
for the Breit-Wigner of the T, 



Pr(s) 



[s - Ml) + zM T r T 



and F~c(q 2 ) is a form factor defined as, 



<T(g, W ) 1^(0)7^(0)10) = F T (q 2 )e*(uj,q) 



(9) 



(10) 



with s*(u, q)^, the polarization four- vector. This form factor can be related 
to the amplitude for T to e + e~ , 



1 2 {^af 
6tt^ 



2 M T 



(11) 



Notice that all the hadronic uncertainties in our process will be included in 
this unique form factor. In Eq. (H) we can see the coupling g\ is only in the 
T^^Za) and T^^Zaj) amplitudes, so, as these amplitudes will contribute 
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to dominant order to the observables through interferences with 7^(7), this 
means that only the P-odd observables will contain the information about 
the g v coupling linearly. Then we are going to analyze the polarizations and 
the P-odd correlations. 

The polarizations of r~ are given as follows 

= P(+,+),(+,+) + P(+ -),(+,-) ~ P(- +),(-+) - P{- ,-),(--) 
= 2 J R e {T (+) _ )il (7)T ( * +i _ )il ( 7 Z A )}(l + cos 2 ^) 
+ 4i?e{T (+ ,_ )il (7)r ( * +i „ )il (Z A7 )}cos0 (12) 



j^P* (0) = P(+,+),(-,+) + + + P(-,-),(+,-) 

= -2v / 2[ J Re{T {+i+)il (7)T ( * +i _ )il ( 7 Z A )}sin^ costf 

+ (Re{T {+>+)>1 (j)T{ +j _ )A (Z A y)} 

+ Jte{T (+f _ )|1 ( 7 )27 +i+)il (Z i i7)})Bintfl (13) 



-*(/>(+,+).(-+) + P(+ -),(-,-) - P(- +),(+,+) - P(- -),(+-)) 

2 v/2 [/m{T (+i+) ,i ( 7 )T ( * +i _ )>1 (Z AjA ) } sin 
2/m{T (+)+))1 (7)T ( * +) _ ))1 (7) } sin 6 cos 6] (14) 

As we can see in these expressions, only P z i and P x / contain information on 
g v , because both are P-odd, T-even observables, but P y > has no information, 
because it is P-even, T-odd. The T-odd observable P y / needs the interference 
between resonant and non-resonant amplitudes, that on the T peak are rel- 
atively imaginary. By the same argument, there is no interference between 
resonant and non- resonant pieces in the T-even observables, like P z > and P x >, 
on the T peak. If we compare the longitudinal polarization, Eq. fll2|), with 
the result we obtained in Eq. (P the difference is in the non-resonant terms 
proportional to gv in Eq. (^). Then these new terms are suppressed on the 
T peak by a factor g v a 2 Qb/ (9g v b.r.(Y — > e + e~) 2 ) ~ 3.5 x 10~ 4 , and so our 
estimate in Eq. ([!]) is very good. This means basically, that this observable 



+ 
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will be, in the forward direction, P z >(9 = 0) ~ .185 ■ g\ ~ 0.032. In Fig. 2, 
we can see a plot of this polarization as a function of 9. 

On the other hand P x > contains similar reduced helicity amplitudes as P z i, 
the only difference is that, as usual, transverse polarizations are suppressed 
by a mass insertion, that is a factor M T /p° ~ .38. Basically, apart from a 
different angular dependence, this is the reason that makes this observable 
less sensitive to g v as we can see in Fig. 3. We get for instance P x >(9 = 
7r/2) ~ —.063 • g v , slightly worse than P z >. 

The information contained in the r + polarizations is closely related to 
that of the r~ 

p(+) _ _ P H 

r z' ~ r z' 

p(+) _ p(-) 

^ +) = (15) 

Finally, we can also consider the spin correlations between r + and r - . 
We define these spin correlations as follows, 



^C tj (9) = Tr(at ] a^ P n (16) 



With this definition is easy to see that the only P-odd observables will be 
the C XtV , C z ^ y , Cy tX and C VtZ correlat ions , 

do , . . . 

C zy {9) = H-p ( + , + ),( + ,-) + P(- +),(- ,-) + P{+ ,-),(+,+) - P(- -),(-,+)) 



= -2v / 2[/m{T (+i+)il (7)T ( * +i _ )il ( 7 Z A )}sin^ costf 
+ (7m{T {+i+))1 ( 7 )r ( * +) _ );1 (Z A7 )} 

/m{r (+ ,_ )>1 ( 7 )7f +>+)>1 (Z i4 7)})8iiie] (17) 

<W0) = (18) 



^Gn/(0) = i(-p( + , + ),(-,-) - P(- +),(+ -) + P(+ -),(-+) + P(- ,-),(+,+)) 

= 27m{T (+i _ )il ( 7 )T ( ; i _ )il ( 7 Z A )}sin 2 e (19) 
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Cyx{9) — —C xy {0) 



(20) 



As we already pointed out and we can see explicitly in Eq. (|I2|)-([2"0|), all 
the relevant observables are associated to the amplitudes T^^jZa) and 
T^^(Za'j), and the relative strength of these observables to the dominant 
term is set by the factor 8GfM^/ '(47tq;a/2) ~ 0.064. Unfortunately, the P- 
odd correlations are also T-odd which means they need an imaginary part. 
With our set of reduced helicity amplitudes, it becomes necessary to have an 
interference between a resonant and a non-resonant diagram to get an imagi- 
nary part. Therefore, these contributions have the same suppression that di- 
agram (1) with respect to diagram (2), that is, a/(3 b.r.(T — > e + e~)) ~ 1/10. 
Notice also that the helicity structure of C zy is very similar to P x i, and so 
it will have the same suppression factors, unlike C xy which is not helicity 
suppressed. 

Then it is clear that our main observables will be the longitudinal polar- 
izations of both r's, then P x > and finally C xy and C zy ordered from the most 
relevant to the least one. 

In the next section we are going to connect these observables to measur- 
able quantities, analyzing the angular distributions of the decay products of 
the two r. 

3 Decay of a polarized r 

The main r decay channels are presented in Table 2. The purely leptonic 
decays have branching ratios of 17.65% for muons and 18.01% for electrons. 
Unfortunately, these decay modes have two neutrinos in the final state, which 
implies we can not reconstruct the r direction. Then, their sensitivity to the 
r polarization is small M compared with the hadronic decays. 

So we will concentrate on the hadronic r decays which have only one 
neutrino in the final state, which allows to reconstruct the r direction if both 
t decay hadronically ||. These decays are r~ — > n~v T , with a branching 
ratio of 11.31%, t~ — > p~v T which corresponds almost exactly to the two 
pions channel, with branching ratio 25.24%, and t~ — > a{v T which is given 
by the sum of the three pion final states. In this work we will concentrate, 
as an example, in the decays r~ — > n~v T and r~ — > p _ z/ r , other r decay 
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channels have been studied elsewhere, for instance r — > a\v T can be found in 
[|T(i|l for LEP physics, 

3.1 T~ — >• 7T~Z/ T 

This channel has been used for a long time to measure the r polarization 
because of its good sensitivity. In this decay, we can easily get the differential 
decay width as, 

f§i = rJ 1+p ' k ^ < 21 > 

where k is the unit vector in the direction of the pion. 

In Eq. (|21p we can see that the polarization effects are not suppressed in 
the angular distribution of the decay pions. This is due to the fact that in 
this process there is only one reduced helicity amplitude, and so the complete 
angular distribution is necessarily proportional to this unique reduced am- 
plitude. So this channel is specially sensitive to the r polarization, and this 
has been the reason for its popularity as polarization analyzer in r decays. 

With these elements, following App. A, we can build a complete two 



steps angular distribution j|, [L1J] for the whole process, e e + — > r r H 

(-K-V T ){-K + V T ) 



dtt dtt + dtt_ dtt 
—C zz cos 6L cos 9 + — C xx sin 6L cos (f)^ sin 9 + cos (j> + 
—C yy sin 6L sin 0_ sin ^ + sin <p + — C za; cos 9- sin # + cos 0+ 
— C xz sin 6'„ cos 0_ cos 6 I + — C ZJ/ cos 6L sin 6 I + sin <f> + 
—C yz sin 6*_ sin 0_ cos ^ + — C^.^ sin 6_ cos 0_ sin 9+ sin + 
— sin 9_ sin 0_ sin 9 + cos0 + ) (22) 

where (6 1 -!-, <fi±) is the direction of the 7r =1= in the rest frame of the and /c-t 
is the unit vector in this direction. 

This will be the cross section we have to study to extract gy from this 
channel. In section 4 we will analyze the sensitivity of the different observ- 
ables we can construct from this cross section. 
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3.2 T —> p V T — > 7T 7T Z/ r 

In the channel r — > pz/ r we have a spin 1 particle in the final state. This 
implies that we have two different helicity amplitudes in this decay, and so 
different combinations of these amplitudes enter in the polarized and unpo- 
larized pieces, 



dV T — ip v 1 



+ pv dCl An 



[1 + a p Pk(n)] (23) 



where a p is a ratio of reduced helicity amplitudes, which we can get in terms 
of the masses as, 



a P \rr 12 , It 1 12 H/T2 , ow2 u -^ du 



KO-1/2 


2 _ 


^-1,-1/2 


2 M 2 - 2M 2 


1^0,-1/2 


2 + 


^-1,-1/2 


2 M 2 + 2M 2 



Then, we can see from here that in spite of its bigger statistics the sen- 
sitivity at this level is smaller than in r — > 7rz/ T . However, as it was pointed 
out in ||, this situation can be improved if, in addition, we try to get some 
extra information on the p helicity To do this, you have to include another 
step in this chain of decays, and analyze the decay p~ — > 7r~7r°. The cross 
section for the whole process, e~e + — > t~t + — > ((7r _ 7r )z/ r )(p + z/ T ), can be 
written as, 



da da , _ , _ 1 dT T - T +^, p - up + D 
- — tee — > r r • : — — - 



dttdtt±dttldtt 2 dtt y ' T 2 dfi^dftj- 



(25) 



where the r~r + decay amplitude to p~vp + u is given below in Eq. (0), and 
we have added the last term which is the decay width of a polarized and 
aligned p into two pions. The expression for this decay width is 



p — t7r7T 



l - [1 - v^£PS(0 2 A,O)i2,iv] (26) 



In this equation we can see that the p polarizations do not appear in the 
decay angular distribution because this is a strong decay and therefore P- 
conserving. The alignments, higher order multipole parameters, appear but 
the polarizations do not. 
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Now, we have to calculate the density matrix for a p coming from the de- 
cay of a polarized r in its center of mass frame, and then apply the necessary 
Wigner rotation [[HJ, as explained in App. A, to get the density matrix in a 
frame where the r is moving. 

This density matrix of a single p from the decay of one of the two rs also 
contains information on the decay of the other r if we study the correlations 
and do not integrate the second r decay. So, as we are interested in the 
measure of correlations between the two rs, we will study the p~ density ma- 
trix in a decay t~t + — > (p~v T )(p + i? T ) and in a decay t~t + — > (p _ i/ r )(7r + z/ r ). 
Naturally this two different density matrices will coincide when we integrate 
completely the r + decay products. 

Following Appendix A, we can write the complete p~ density matrix from 
a decay t~t + — > (p~ v T )(p + u T ) as 

Pix-fj,'- = Sa + A„A^_A'_ Yjv-v'_v + /(t^l/2)A+(^+)^M-' / -( a; -)/(i_-l/2)A_(^-) 

P\X-,\ + ),{X'_,X' + )f[v + l/2)\' + (^+) rf M>l - 1/2)A'_ (^-) ( 27 ) 

where we have used a reference system in LAB with the r~ in the z-axis 
and the initial beams in the x-z plane, to simplify the expressions for the 
Wigner rotations. 

We define an effective density matrix, p without the Wigner rotations 
that, if we integrated completely the fl + variables, would correspond to the 
density matrix in the r~ rest frame, which is 

/W_ = d M-v-( u -)Pv-vLdv,'_v'_{w-) (28) 

The next step is to calculate this effective density matrix in terms of 
reduced helicity amplitudes. In a decay we define the helicity amplitudes as, 



MO, <P) = f-^V{* Ui _ V2 {^ 6, 0)T£ (29) 



With this definition and Eqs. Q2"7D and (28), we get our effective density 
matrix. In this process, r~ — > p~v T , we have only two reduced helicity 
amplitudes if we take M v = 



KT_ lt . 1/2 = %AG F V ud F p (q 2 ) yj M T k 2 



KT ^ 1/2 = t2V2G F V ud F p (q 2 ) {mJ 2 ^ (30) 
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and also two amplitudes in the decay r + — > p + u T1 

KT hl/2 = -iAG F V ud F p (q 2 ) JmJ 2 

KT , 1/2 = t2V2G F V ud F p (q 2 ) JmJ 2 ^ (31) 

where F p (q 2 ) is a form factor defined as, 

(p(k,u)\T\n(k 2 )n(-k 2 )) = F p (q 2 )e*(u,q) p k^ (32) 
The rest of the notation is self-explanatory. The total amplitude p — > tttt is, 

Tp^r = ^~d T -i -1/2I 2 + l^o,— 1/2 1 2 ) (33) 

Following App. A, we now calculate the CM. multipole parameters cor- 
responding to the density matrix p in terms of these reduced amplitudes, 
taking into account that we do not integrate the direction of the second p. 



-^d f = Tr{p} = p^ 1 + p 0fi 

= (\T-1 ,—1/2 1 2 + I ^0,-1 /2 1 2 ) ( I ^1,1/2 1 2 + 1^0,1/2 1 2 ) 

[1 + a p P ■ k~ — aP ■ k + — a p a(C zz cos cos 6 + 
+ C xx sin 6_ cos 0_ sin 6 + cos <p + + C yy sin 6_ sin 0_ sin 6 + sin <f> + 
+C ZX cos 6L sin 9 + cos (f) + + C xz sin 6L cos cos 9 + 
+ C zy cos 6L sin 6* + sin (f> + + C yz sin 6L sin 0_ cos 6 I + (34) 
+ C XJ/ sin 6L cos 0_ sin # + sin + + C yx sin 6L sin 0_ sin 6> + cos + )] 

here (9±, <p±) is the direction of the p ± in the rest frame of the and k± is 
the unit vector in this direction. a p was defined in Eq. (S) and a is equal 
to a p if the r + decays to p + v T1 but is equal to 1 if it decays to ix + u T . This 
definitions hold for all the observables that follow. 



Tr{p}i 20 = V 15 " 2 /M = V^(l T -i--i/2| 2 + l^o,-i/ 2 | 2 ) 
(|T 1>1/2 | 2 + |r , 1/2 | 2 )[ 7p - f3 p P- ■ k- - a lp P+ ■ k+ 
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+af3 p (C zz cos 9- cos#+ + C xx sin 9- cos 0_ sin # + cos <j) + 
+C yy sin 6L sin 0_ sin # + sin + + cos sin 9 + cos + 
+Ca; 2 sin 9- cos 0_ cos 9 + + C 2J/ cos sin 9 + sin + 
sin 9- sin 0_ cos 9 + + C XJ/ sin cos 0_ sin # + sin <p + 

+C yx sm9_ sin0_ sin6> + cos0 + )] (35) 

where we have introduced two new coefficients f3 p and 7 P defined as 

|T_ 1 _ 1 / 2 | 2 -2|T _ 1 / 2 | 2 2M 2 -2M 2 / x 

7p= W 2 7 i J 7 = „2 , ^ =-1-18 (36) 



l^0,-l/2| 2 +| 
l^-l,-l/2| 2 + 


^ —1,-1/2 1 

■ 2|T 0; _i/i 


! M 2 + 2M 2 
,| 2 2M 2 + 2M 2 


KO-1/2 


| 2 + 


l T -i,-i/2 


| 2 M 2 + 2M 2 



| T 1/2I 2 -' ! ' ,: ' " ' r ' 1 r> 

we also need t 2 i and t 22 , 



Tr{p}i 2l = y/Zp-w = Jf Q (\T_ lt _ 1/2 \ 2 + |T ,_ 1/2 | 2 ) 



— i sin 0_) 
#4 



J 7 t^21 = Y 10^-1,0 = V 10 U J -1,-1/21 ^1^0,-1/21 ; 

(|Ti,i/ 2 | 2 + |T 0il/2 | 2 )[5 p (P7(-sin^) + p-(cos0_ cos0_ - i 
+P~(cos0_ sin0_ + % cos0_)) — 5 p a(C zz (— sin#_) cos^ + 
+ C xx (cos 9- cos <fi- — i sin 0_) sin # + cos + 
+C yy (cos sin 0_ + % cos 0_ ) sin 9 + sin + + C zx ( — sin ) sin 9 + cos 0. 
+C xz (cos9- cos0_ — isin0_) cos6> + + C zy (— sin#_) sin6> + sin0 + 
(cos 6L sin 0_ + % cos 0_ ) cos 9 + 
+ (cos cos 0_ — i sin 0_ ) sin 6'+ sin <p + 
+C yx (cos sin 0_ + % cos 0_) sin 9 + cos 0+)] 



Tr{p}t 22 = y^P-M = ( 39 ) 



we have also introduced a new coefficient 

J = r -^-V^.-V2 = ^MeMi = 445 (40) 
P |T ,_ 1/2 | 2 + |T_ W/2 | 2 M 2 + 2M 2 U ' 445 



In general, these multipole parameters can be complex, as t 2 ±, except the 
parameters with L = 0, which are always real. 
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It is very important to notice that unlike a p and 8 P , that are small, both 
7 P and j3 p are bigger than one and so they can enhance the information on 
the r polarizations. 

As we have already pointed out, the only difference between the decays 
t~t + — > (p~ v T )(p + v T ) and r~r + — > (p~z/ r )(7r + z/ T ) at this level is the value of 
the coefficient a 



t + — > p + u T => a = .456 

r+ -> 7T+z/ T =>a=l (41) 



The final step to get the alignments appearing in Eq. (j26|) and Eq. (P5|) 
is to make the Wigner rotation, Eq. ( |A.13| ). 



f-20 



L 2M 



l 20\ 2 



COS CJ_ — 



sin cos 



(42) 



' 21 



C ^1J\/( U; -)^2M 



-£20 cos w - s i n w - + ilm{t 2 i } cos u;_ 
+_Re{t 2 i}(cos 2 o;_ — sin 2 c<j_) 



(43) 



^22 = d\ M (uJ)t 2 M = ho^f sin 2 cj_ - ilm{t 2 i} sinu- 

—Re{i 2 i } cos uj_ sin cj_ (44) 

where us is the Wigner rotation associated with the boost from the r rest 
frame to the e + e~ CM. frame which transforms the p four-momentum, k' p = 
Ak p . These rotations have the following expression 

M p \p T \sm6 E p E' p M T -ZM 2 p 

smu = — - , coscj = =r- — =; (45) 

M T \k' p \ | k p || k' p | M T 

with (Z,p T ) the r four-momentum in the e + e~ CM. frame and 9 the angle 
between the p and the direction of the boost in the r rest frame. 

Now we have the complete angular distributions of the different chains 
with final states: 
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• (tt V T ){lt+V T ) 

• (71^71°u t )(ti + u t ) + (7r + 7r z/ T )(7r~z/ r ) 

• (7r _ 7r°z/ T )(p + z/ r ) + (7r + 7r°z/ T )(p~z/ r ) 

Next we are going to find what is the statistical accuracy one can obtain 
in the measure of g\ 



4 Statistical sensitivity 

To estimate the obtainable precision in the measurement of g\ in a B-meson 
facility, we use the formalism of references [T3|, In these references 



they call "ideal statistical error" of a parameter p which enters a function 
f(x, y;p) to be determined experimentally, to the error obtained from a least- 
squares-fit to this function with N events. To obtain this error, we use that 
for large number of events, N, the likelihood function approaches a gaussian. 
Then, if the function f(x,y;p) is normalized to one on the physical region, 
the ideal statistical error is given by 

°* = i^/ ^ ~^~^p" ' ,Xn,P h 2 ■ f( x u ■~,x n ;p)dx 1 . . .dxn}- 1 (46) 

The word "ideal" stands for the fact that we are not considering the 
efficiency of the detectors, effects of finite experimental resolution and we 
assume an ideal distribution of the N events according to f(x,y;p). 

In this case our function f(x,y;p) will be the normalized cross section 
to the different channels. First we study the reaction e~e + — > t~t + — > 
(7r~ u T ) (tt + V t ) . Our result is, in this channel, 

cv = —,= (47) 



where N is the number of e~e + — > r _ r + — > (tt~ v T ) (tt + v T ) events. This means 
that in a B-meson facility with 10 8 T(1S) produced per year one could get 
a sensitivity to g\ of 6 • 10~ 2 only with this channel. In Eq. (^) is all 
the information available in the process, but not all these observables will 
be useful for our purposes. In particular the P-even, T-even correlations 
(C zz ,C xx ,C yy ,C zx ) get contributions from |T^(7)| 2 , so they are order 1 and 



15 



they have no information on g v . On the other hand, the P-odd, T-odd 
correlations are not as sensitive to this parameter as the polarizations, then 
we can ask whether we can improve the sensitivity by integrating out some 



of the final state variables. From Eq. (|22|) we would like to eliminate the 
P-even correlations maintaining the polarizations and P-odd correlations, 
but unfortunately we can not achieve this result integrating out some final 
state variables. The only interesting possibility is to consider both r decays 
independently, integrating out one of the two pion directions with respect 
to the t. In this case we are considering the e~e + — > t~t + — > ({ii~v t )t + + 
t~(tt + v t )) events, where both r are hadronically reconstructed. By doing this 
we increase considerably the number of events, because now we include also 
the events in wich the second r decays to p and a%. This is approximately a 
factor of 5 in the number of events for each tau, but we have two independent 
events for each T decay, then the number of events increases a factor of 10. 
The sensitivity one can get with this new decay distribution ( we obtain it 

it (e+,(f> + ) ) is 



directly from (f22|) integrating out (#+,</>+) ) is 

14.6 



and now N is the number of e~e + — > t~t + — > ((h~v t )t + + t~(it + i> t )) events. 
We can see that the difference between Eqs. (fTTD and P5| ) is roughly a factor 
v2, this is due to the fact that in Eq. ( ]47j ) we included both the r + and t~ 
polarizations while in Eq. ( |4~8|) we take into account only one of them. On 
the other hand, this means that the P-odd correlations do not improve the 
measurement of g v in a significant way. Again with 10 8 T per year, one can 
get a sensitivity of 2.3 • 1CT 2 . Even more, the simplest polarization analyzer 
we can use is the energy of the pions. The energy of the pions in LAB is 
related with the angle in CM. of the r 

K - E ^ + ^ 9 - (49, 

where (E*, A;*(0_)) is the four-momentum of the pion in the r CM., and 
(E T ,q) the four-momentum of the r in LAB. Again from Eq. (0), if we 
integrate all the angular variables but 9- and make this change of variable 
to E v one gets a sensitivity of 
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9 - s/N 

And, in this way, we do not put any restriction on the second r decay, then 
with 10 8 T per year one can get a sensitivity to g v of 2.7 • 10~ 2 , simply using 
the 7rv channel and measuring only the pion energy. From this point of view, 
it is evident that, if we consider only the irv decay channel, this is the best 
strategy to measure g v , because one can use all the r — > irv events and is 
experimentally simpler. 

We have also studied the channel r — > pv as a polarization analyzer. 
Then if we apply Eq. fl4"6D to the complete distribution e~e + — ► t~t + — > 
((n~TT°)v T )(p + v T ) + ((tv + tv°)v t )(p~v t ) given by Eq. (|25|) one gets 

12.4 . . 

a a b = —= (51) 
9v s/N 



that is slightly worse than our result in the e~e + — > t~t + — > (tt~ v T ) (ir + v T ) 
channel, but now we have more events because these channel has a bigger 
branching ratio. This translates in a B-meson facility with 10 8 T produced per 
year in a sensitivity to g v of 2.1 ■ 1CT 2 . Other possibility is to use a combined 
channel as e~e + — > t~t + — > (((7r _ 7r )z/ T )(7r + P T ) + (7r~z/ T )((7r + 7r )z/ T )) where 
the sensitivity to the polarization of the r + — > tt + v t is better than r + — ► p + u T 
if we do not analyze the next decay p + — > 7r + 7r . The result is, 

9.1 , . 

a qb = —= 52 

and the number of events is similar to the previous case. The sensitivity to 
g v that one can reach in this channel, with 10 8 T, is 2.3 • 1CT 2 . As in the 
7r channel we can increase the statistics by integrating one of the decays, 
although we lose some sensitivity. We integrate the direction of one of the p 
or equivalently the ir direction in Eq. (p5|) and only require this r to decay 
hadronically. Then we keep simply the decay of a r to pv and then the decay 
of p to 7T7T. Then we get a sensitivity of, 

°9\ = -7^ (53) 



but now the number of events has increased a factor of 10, because N is the 
number of e~e + — > t~t + — > (((7r _ 7r°)z/ T )r + + t~ ((7t + tt )v t )) events. Again 
with 10 8 T per year, one gets a sensitivity of 1.6 • 10 -2 . 
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After this analysis we can combine a series of independent measurements 
into a final value for g v with an error given by 

* = (E -V 1/2 (54) 

i a i 

then we can combine the error obtained with the channels e~e + — > 
(7r~z/ T )(7r + z/ T ), e~e + — > (((7r~7r°)z/ T )(p + z/ r ) + ((p~v> T ) ((7r + 7r°)z/ T )) and e~e + — > 
(((7T _ 7r°)z/ r )(7r + z/ T ) + (7T _ z/ T )((7r + 7r°)z/ T )) and with 10 8 T one gets a sensitivity 
to g h v of 1.5 • 1(T 2 . 

On the other hand, we can also combine the errors obtained in the e~e + — > 
(({it~ -k )u t )t + +t~ ({it + it°)u T )) ande _ e + — > ((7r~z/ T )T + +T~(7r + z/ T )) channels, 
and again with 10 8 T one gets a sensitivity of 1.3 • 10~ 2 . 

Notice this result has been obtained with a sample of 10 8 T, for a different 
number of T produced, the sensitivity to g\ would simply re-scale by a factor 
J 10 s /N r . 

5 Conclusions 

In this work we have studied the possibilities of a high luminosity B-meson 
facility to measure with high precision the Z — bb vector coupling. At the 
energies of T(IS') we have used the r~r + channel to determine this coupling 
through the r polarizations. A complete analysis of the hadronic decay modes 
of the r lepton has been done, with special attention to the r~ — > 7r~z/ T and 
t~ — > p~v T as polarization analyzers. We have built the complete correlated 
cross section with the decays of both rs and from here we have calculated 
the ideal statistical errors obtainable in the measure of g v . We have found 
that in a one year run in a B-meson facility with 10 8 T per year, one can 
get a sensitivity of 1.3 • 1CT 2 , comparable with the present precision in this 
coupling from the LEP/SLC measurements of Rb and A b . 
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APPENDIX A 

In a decay, A(a) — > B(Xb) + C(A C ), we can obtain the density matrix, p out , 
describing the complete final state of particles B and C in terms of helicity 
amplitudes and the initial density matrix in the following way, 

pfh = E k,„ («0 *&(n) ft (j> (no (a.i) 



where A = (A&, A c ). 

It is very convenient to express our initial density matrix in a basis of 
irreducible tensor operators, Tl,m, with coefficients tL,Af, @, 01, 



2/ 



rm^WTTU 2L + 1)t ^ (A ' 2) 

these coefficients, t^ Ml are the so-called multipole parameters, 

Tr(p in )t^ M (9) = Tr{p™T^ M ) = J2(p° u %yC(lLl\aMa>) (A.3) 

a, a' 

where these C(jLj\m'Mm) are Clebsch-Gordan coefficients. 

For L = 1, we can relate, @, the usual polarizations and the multipole 
parameters, 

P*> = - Vi) = "2 Re[h t i] (A.4) 

Py = + *i,-i) = -2 Jm[ti,i] (A.5) 

P a i = v/2 £ lj0 (A.6) 

Using Eq. ( |A.2|) to express the initial density matrix in the basis of 
irreducible tensors and replacing the helicity amplitudes in the CM. frame 
of the decaying particle with Eq. ( p9|) , we get the final density matrix in 
terms of the multipole parameters of the A particle and the reduced helicity 
amplitudes, 



p$ = Trin^^Y. v / 2TTl (-y- X 'T x Tl (A.7) 



L,M 

41/ C(j,j, L\X', -A, A' - A) V$*(Q) 
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This density matrix contains all the information available in the process. For 
instance, if we want the angular distribution we just have to take the trace 
on A, and in the same way the density matrix for one of the final particles is 
obtained taking the trace on the helicities of the other particle. 
To define a complete set of observables we generalize Eq. ( |A.2|) 



P 



out 



Tr(p' 



,out\ 



(2j 1 + l)(2j 2 + i; 



(b) rp(c) 



M' 



2j 

£ (2L+l)(2L' + l)C2 >M)i , M ,T J 

L,M,L',M' 

(A.8) 
include all the 



and then these generalized multipole parameters, Cl,m,l',m', 
information on the density matrices of particles B and C, t^ M 



C 



L,M,0,0; 



t 



L M — Co,o,l,m and additionally the correlations between them. 

It is very important to notice here that Eq. ( |A.7p is only valid in the 
CM. frame of the decaying particle. However, in general in the LAB frame 
the decaying particle will move with a momentum different from zero. So, we 
will be interested in the transformation properties of these density matrices 
under Lorentz boosts. 

The transformation of an helicity amplitude under a boost from the CM. 
frame to the LAB frame is just a rotation, the so-called Wigner rotation 

= 0, is [HI, 



l|, [Tg, [lTfl , that if we choose 



(CM) 



f\ h \, 



,0 = 0) 



(-) 



A'-A c 



Ale 

A' Ar 



where this rotation, of angle Ui, is given by 0, 
Mi sinh k sin 9 



sin LOi 



COS LOi 



Pi 



EtEi - cosh kM? 
I Pi 1 1 Pi I 



= 0) 
(A.9) 

(A.10) 



where k is the parameter of the boost from the CM. frame of the decaying 
particle to the LAB, related to the velocity by v — tanh k. Mi the mass of the 
particle and (Ei,pi(9)) its four-momentum in CM. of the decaying particle, 
transformed under the boost as p[ = Api, in a frame where the boost is along 



the z-axis. From Eq. (|A.10| ) we can see that under a boost collinear to the 
particle three-momentum, 8 = 0, our states do not suffer any rotation. 

The helicity amplitudes in Eq. (|A.9| ) are functions of two variables. For 
instance, we could choose the invariant variables (s,t,u), keeping the same 
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expression in any frame. Nevertheless, as we have seen in Eq. (|29|), these 
helicity amplitudes have a specially simple form in terms of CM. variables. 
Then, we have used this freedom to express, both the CM. and the LAB 
helicity amplitudes in Eq. (|A.9|) in terms of CM. variables. 

As can see in Eq. ( |A.1|) , density matrices are a product of two helicity 
amplitudes, this means that using Eq. ( |A.9|) we can get the transformation 
of the density matrix of particle B 



p™ = d»{ Ub )-p tJM (A.ll) 

It is very interesting to see how the multipole parameters are affected by 
these rotations. We use Eq. ( |A.3| ) to express p CM in terms of these multipole 
parameters 



P L J B = 2J+1 ^lm(2L + 1) Eav ^MijLj^'MX^/Mt^; 

= ^1 J2lm(2L + l){jLj\a>o - a')df}J M {u)t%* (A.12) 

and comparing again with Eq. (|A.3| ) that in LAB we get a new set of multi- 
pole parameters that are obtained simply applying the rotation to the CM. 
ones. 

^lm = y^;d MM ,(u)b)t LM > (A. 13) 

M' 

This is all we need to obtain the multipole parameters and density matrices 
in the LAB frame. 



APPENDIX B 

In this appendix, we present the general method to calculate reduced helicity 
amplitudes, and we apply it to some examples in the e + e~ — > t + t~ processes. 

Reduced helicity amplitudes are easily calculable by means of Eq. (Q) 
from the helicity amplitudes. So, our first step will be to obtain these he- 
licity amplitudes from the Feynman amplitudes we can calculate from the 
diagrams with Feynman rules. Taking into account the normalization we have 
defined for our reduced helicity amplitudes in Eqs. (§) and (|29|) the differ- 
ence between them and the Feynman amplitudes will just be a g 2 -dependent 
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phase space factor, irrelevant in all our observables. So, we can simply define: 

M aM , x _(9) = Kf aX (9) (B.l) 

with M Ct \ +j \_{9) the Feynman amplitudes. 

Now we will explicitly apply this procedure to the calculation of the re- 
duced helicity amplitudes 7(a,a'),£(7) and T^x^^jZ^) corresponding to di- 
agrams (1) + (2) and (3) + (5) in Figure 1. 

The Kinematics in the CM. frame of the e + e~ system is defined by 

It = (£,0,0, |/t) q» = (L+l + r (B.2) 

11= (E,0,0,-\l\) k" = (k°, |^| sin 0,0, \k\cos6) 

where l± is the four-momentum of the e ± , whose helicities are £± = ±1/2, 
the four-momentum of the r~ and the helicities of the r ± will be denoted 
as A± = ±1/2. 

The Feynman amplitudes corresponding to these diagrams are 

mLmm-W = 4^ + 4^l F T(e a )r^(e 2 ))K'(A-,A + ,%^(e-^ + ) 

(B.3) 

M l Z M^°) = ^9A(jQ b g b v \Fr(q 2 )\PAq 2 ) ~ 9v) 

Al{\^\ +I e)g^{^i + ) (B.4) 

where we have followed the notation of sections 2, and we have introduced 
the matrix elements of the leptonic currents, 

Vf(A_,A + ) = u(p_,A_) 7 Mp+,A + ) 
Af (A_, A+) = ft(p_, A_)7"T5V(P+, A+) (B.5) 

Now we need to obtain an explicit expression for these matrix elements. 
To do this, we follow the method of reference |Jl8| , which permits the cal- 
culation of these amplitudes using standard trace techniques. Then, the 
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complete results for the vector and axial currents, in the CM frame and with 
the momenta along the z-axis are, 



VT(A_,A + = A_) = (0,0,0,-2MO 
V^(A_,A+ = -A_) = (0,AE P \-,-2E p i,0) 
Af(A_, A + = A_) = (-4MA-, 0, 0, 0) 
Af (A_, A+ = -A_) = (0, 2p, -4zpA_, 0) (B.6) 

With these results we have all the leptonic currents we need, because they 
are perfectly behaved Lorentz vectors or axial-vectors. Then, we just have 
to rotate them, if the momentum is in a different direction. 

With all these elements, we simply use Eq. ( |B.1| ) and Eq. (Q) with our 
Feynman amplitudes, Eqs. ( P-3| ) and ( |B.4| ), to obtain the reduced helicity 
amplitudes. For instance, the expression for M\ - /o \(0) 

A — (1/2, 1/2), £ — (1/2, — 1/2) 

2 2 

^(V2, 1 /2),(i/2,-i/2 ) (^) = ijCi + 7Q^T(g 2 )| 2 Px(g 2 )) 

(0,-2M T sin^,0,-2M r cos^) • (0, -2^, -2i^, 0) r = 

4(! + 7^l^(g 2 )| 2 Px(g 2 ))4M T ^ (-V2) (^) (B.7) 

where we have applied a rotation to the leptonic current of the r with respect 
to Eq. ( |B.6| ) and we have taken the complex conjugate of Eq. ( |B.6| ) to obtain 
the electron current. The extra minus sign is due to the metric g^ v . In Eq. 
( |B.7| ) we just have to remove the rotation matrix element d w (6), which is 
exactly the last term in this equation. This procedure has to be repeated 
with all the amplitudes and then, finally we get the following results 



KT {+>+);1 (j) = -iW2 - (1 + ^Ql\F r \"P T ) M T ^ (B.8) 

S S — 

^ (+r)il ( 7 ) = -i8 - (1 + -Qi\F r \ 2 P r ) p° ^ (B.9) 
s s 2 

KT {+i+)A (<yZ A )=0 (B.10) 
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KT^^Za) = -i8 ^ g A (^Q b g b v \F r \ 2 P r - g v ) \p\ ^ (B.ll) 
In the same way, we can obtain all the reduced helicity amplitudes in Eq. 
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Figure Captions 



Figure 1: Resonant and non-resonant diagrams for the process e + e 
t + t~ at the T region. 

Figure 2: Longitudinal r~ polarization, P z >{9) 
Figure 3: Transverse r~ polarization, P x >(9) 



27 



Table Captions 



Table 1: Dominant T(IS') decay channels. 
Table 2: Dominant r decay channels. 
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Table 1 



T(1S) I G (J FC ) =0-(l— ) 


Mass My = 9460.37 ± 0.21 MeV 
Width T = 52.5 ± 1.8 KeV 


Decay modes 


Fraction Ti/T 


T + T - 
e + e~ 

J/ip{lS) anything 
7 2h+2h~ 
7 3h + 3h' 
7 Ah+Ah- 


(2.97 ±0.35)% 
(2.52 ±0.17)% 
(2.48 ±0.07)% 
(1.1 ±0.4) x 10~ 3 
(7.0 ± 1.5) x 10~ 4 
(5.4 ±2.0) x 10~ 4 
(7.4 ±3.5) x 10~ 4 



Table 2 



T J=\ 


Mass M T = 1777.001^7 MeV 
Mean life r = (291.0 ± 1.5) x 10~ 15 s 


Decay modes 


Fraction Ti/T 


e~v e v T 
ix~v T 

h~2n°v T 
h~h~h + v T 


(17.35 ±0.10)% 
(17.83 ±0.08)% 
(11.31 ±0.15)% 
(25.24 ±0.16)% 
(9.50 ±0.14)% 
(9.80 ±0.10)% 
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Fig. 1 



